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ABSTRACT 

Context. IRAS 20343+4129 was suggested to harbor one of the most massive and embedded stars in the Cygnus OB2 association, 
IRS 1, which seemed to be associated with a north-south molecular outflow. However, the dust emission peaks do not coincide with 
the position of IRS 1, but lie on either side of another massive Young Stellar Object (YSO), IRS 3, which is associated with centimeter 
emission. 

Aims. The goal of this work is to elucidate the nature of IRS 1 and IRS 3, and study their interactions with the surrounding medium. 
Methods. The Submillimeter Array (SMA) was used to observe with high angular resolution the 1.3 mm continuum and CO (2-1) 
emission of the region, and we compared this millimeter emission with the infrared emission from 2MASS. 

Results. Faint millimeter dust continuum emission was detected toward IRS 1, and we derived an associated gas mass of ~ 0.8 Mq. 
The IRS 1 Spectral Energy Distribution (SED) agrees with IRS 1 being an intermediate-mass Class I source of about 1000 Lq, whose 
circumstellar material is producing the observed large infrared excess. We have discovered a high-velocity CO (2-1) bipolar outflow in 
the east- west direction, which is clearly associated with IRS 1. Its outflow parameters are similar to those of intermediate-mass YSOs. 
Associated with the blue large-scale CO (2-1) outflow lobe, detected with single-dish observations, we only found two elongated low- 
velocity structures on either side of IRS 3. The large-scale outflow lobe is almost completely resolved out by the SMA. Our detected 
low-velocity CO structures are coincident with elongated H2 emission features. The strongest millimeter continuum condensations 
in the region are found on either side of IRS 3, where the infrared emission is extremely weak. The CO and H 2 elongated structures 
follow the border of the millimeter continuum emission that is facing IRS 3. All these results suggest that the dust is associated with 
the walls of an expanding cavity driven by IRS 3, estimated to be a B2 star from both the centimeter and the infrared continuum 
emission. 

Conclusions. IRS 1 seems to be an intermediate-mass Class I YSO driving a molecular outflow in the east-west direction, while 
IRS 3 is most likely a more evolved intermediate/high-mass star that is driving a cavity and accumulating dust in its walls. Within and 
beyond the expanding cavity, the millimeter continuum sources can be sites of future low-mass star formation. 
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1. Introduction 

On the northeastern side of the Cygnus OB2 association, and 
at 1 .4 kpc of distance from the Sun (Le Duigou & Knodlseder 
2002; Sridharan et al. 2002) there is a rimmed feature bright 
at centimeter wavelengths (Carral et al. 1999) and in the H2 
emission line at 2.12 fim (Kumar et al. 2002), which harbors 
at its center the source IRAS 20343+4129. This IRAS source 
is a high-mass protostar candidate of 3200 Lq (Sridharan et al. 
2002) embedded in dense gas (Richards et al.1987; Miralles et 
al. 1994; Fuller et al. 2005; Fontani et al. 2006). When observed 
with high angular resolution, two bright nebulous stars, IRS 1 
(north) and IRS 3 (south), are found inside the IRAS error ellipse 
(Kumar et al. 2002). Either or both sources might account for the 
bulk of the total luminosity of the IRAS source. Comeron et al. 
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(2002) carried out a study of the red and massive objects in the 
entire Cygnus OB2 association and conclude that IRS 1 may be 
one of the most luminous and deeply embedded members of the 
OB association, based on 2MASS JHK photometry and spec- 
troscopy in the range 1.5-2.4 fim. Extending southward from 
IRS 1 and surrounding IRS 3, there is H2 line emission in a fan- 
shaped structure which corresponds very well with a blueshifted 
CO (2-1) lobe detected in single-dish observations by Beuther 
et al. (2002b). The CO structure found by these authors is bipo- 
lar, and the redshifted lobe is centered on IRS 1, suggestive of 
a north-south molecular outflow. Kumar et al. (2002) detect ex- 
tended H2 emission in the east- west direction toward IRS 1, and 
attribute this emission to arise in a circumstellar disk, perpendic- 
ular to the north-south outflow. Although all these observations 
seem to indicate that IRS 1 is a high-mass YSO, no significant 
amount of ionized gas is found associated with IRS 1 (Carral et 
al. 1999). The only compact centimeter continuum source in the 
region is associated with IRS 3, which is interpreted as an Ultra 
Compact H n (UCH 11) region ionized by a B2 star (Miralles et 
al. 1994; Carral et al. 1999). Single-dish observations at 1.3 mm 
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reveal two millimeter continuum peaks lying on either side of 
IRS 3 (Beuther et al. 2002a; Williams et al. 2004). Thus there 
is no clear evidence on whether IRS 1 and/or IRS 3 is the in- 
frared source producing most of the IRAS luminosity, and their 
relation with the dust condensations and the outflow emission 
remains unclear. 

In this paper we present SMA observations of the contin- 
uum emission at 1.3 mm and the CO (2-1) emission toward 
IRAS 20343+4129. Our data provide an angular resolution six- 
teen times better in area than that of the single-dish observations, 
allowing us to gain insight into the nature of each source in the 
region. 



2. Observations 

The SMA!!] (Ho et al. 2004) was used to observe the 1.3 mm 
continuum emission and the CO (2-1) (230.53796 GHz) emis- 
sion toward IRAS 20343+4129. The observations were carried 
out on 2003 August 3, with 6 antennas in the array. We note that 
this was one of the first days for the SMA to work with 6 anten- 
nas. The phase center was ff(J2000) = 20 h 36 m 07^3, <5(J2000) = 
+41°39'57'.'20, and the projected baselines ranged from 13.1 to 
119.8 m. The pads of the antennas were 1, 4, 5, 8, 11, and 16, 
which correspond to a hybrid between the compact and extended 
configurations. System temperatures were around 200 K. The 
full bandwidth for each sideband at that time was 0.984 GHz, 
and each sideband was divided into three blocks with four base- 
bands in each block. The correlator was set to the standard mode, 
which provided a spectral resolution of 0.8125 MHz (or 1.06 
km s" 1 per channel) across the full bandwidth of 1 GHz. The 
FWHM of the primary beam at 230 GHz was -45". 

The raw visibility data were flagged and calibrated with the 
MIR-IDL0 package. The passband response was obtained from 
observations of Uranus, which provided flat baselines when ap- 
plied to Neptune. The baseline-based calibration of the ampli- 
tudes and phases was performed by using the source 2015+371. 
Typical rms of the phases was ~ 60°, yielding a positional un- 
certainty of ~075. Flux calibration was set by using Uranus, and 
the uncertainty in the absolute flux density scale was ~20%. 

Imaging was conducted using the standard procedures in 
AIPS (for the line), and MIRIAD (for the continuum; Sault 
et al. 1995). The channel maps were cleaned adding a value 
for the zerospacing parameter of 100 Jy (estimated from the 
single-dish observations of Beuther et al. 2002b) in the IMAGR 
task of AIPS, and using different clean boxes. The continuum 
map was made using only the lower sideband and excluding 
the CO (2-1) line (the upper sideband was noisier and com- 
bining both sidebands did not result in a higher S/N). We used 
three clean boxes to clean the continuum, one box for each con- 
densation detected with single-dish, and the other box toward 
IRS 1 . Cleaning with these boxes minimized the negative side- 
lobes in the final cleaned map. The final synthesized beam is 
3'.' 47 X 2'.'64, with P.A.= -38°0, the rms noise of the continuum 
map is 2 mJy beam -1 , and the rms noise of the 2.11 km s" 1 wide 
channel maps is 0.4 Jy beam" 1 . 



1 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 

2 The MIR cookbook by C. Qi can be found at 
http://cfa-www.harvard.edu~cqi/mircook.html 



Table 1. Parameters of the sources detected above 5<x at 1.3 mm 
in the IRAS 20343+4129 region 



Position" /r akb S h v Mass c 



Source 


a(J2000) 


(5(J2000) 


(mJy beam 1 ) 


(mJy) 


(Mq) 


MM1 


20:36:05.56 


+41:40:00.2 


26.8 


44.8 


1.2 


MM2 


20:36:06.30 


+41:40:00.7 


40.5 


40.5 


1.0 


MM3 


20:36:06.31 


+41:39:56.4 


29.9 


34.0 


0.9 


MM4 


20:36:06.62 


+41:40:00.6 


39.7 


45.3 


1.2 


MM5 


20:36:07.49 


+41:40:12.8 


25.8 


34.8 


0.9 


MM6 d 


20:36:07.56 


+41:40:08.0 


16.9 


31.5 


0.8 


MM7 


20:36:08.19 


+41:39:54.8 


16.7 


28.0 


0.7 



a Positions corresponding to the intensity peak. 
b Corrected for the primary beam response. 

c Masses derived assuming a dust temperature of 30 K, and a dust 
mass opacity coefficient from Ossenkopf & Henning (1994, see 
main text). The uncertainty in the masses due to the opacity law 
is estimated to be a factor of 4. 

d Associated with IRS 1. 



3. Results 

3.1. Continuum 

Figure Q] shows the continuum emission observed with the SMA 
at 1 .3 mm. The emission is found basically toward three posi- 
tions in the field: to the west of IRS 3 (where we found the 
strongest condensations, detected up to 18<x), to the east, and 
to the north of IRS 3 (where we found condensations up to 9<x). 
The condensations to the west and to the east of IRS 3 are co- 
incident with the single-dish peaks of emission (Beuther et al. 
2002a), and are not associated with any infrared source. The 
emission to the west of IRS 3 contains four condensations, MM1 
to MM4, with a total flux density of 230 mJy (corrected for the 
primary beam response), while the emission to the east has one 
condensation, MM7, of only 28 mJy. This is different from the 
single-dish measurements, for which the eastern condensation is 
stronger than the western condensation by almost a factor of 2. 
We estimated that the SMA has picked up only 2% of the flux 
density of the single-dish eastern condensation (of ~ 30" in size), 
and 30% of the flux density of the single-dish western condensa- 
tion (of ~ 20"). Thus, the emission of the eastern condensation 
is essentially extended, and has been filtered out by the SMA, 
while the emission in the western condensation consists of dif- 
ferent compact millimeter sources. 

In Table Q] we list the position, peak intensity, flux density 
and mass for each millimeter continuum condensation detected 
above 5<x. In the derivation of the masses, we assumed that all 
of the continuum emission is dust emission which is optically 
thin, and adopted a gas-to-dust mass ratio of 100 and a dust 
mass opacity coefficient at 1.3 mm of 0.9 cm 2 g" 1 (agglomer- 
ated grains with thin ice mantles in protostellar cores of densi- 
ties ~ 10 6 cm" 3 ; Ossenkopf & Henning 1994). There is a fac- 
tor of 4 in the uncertainty of the masses due to uncertainties 
in the opacity law. As for the dust temperature, we found two 
estimates in the literature. From NH3 observations toward this 
region, Miralles et al. (1994) derive a rotational temperature of 
~ 20 K, which can be considered a lower limit for the dust tem- 
perature, as gas from protostellar envelopes is mainly heated by 
collisions with warm dust grains (e.g., Ceccarelli et al. 1996). 
On the other hand, by fitting two greybodies to the spectral en- 
ergy distribution of IRAS 20343+4129, Sridharan et al. (2002) 
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Fig. 1. White (positive) and thick grey (negative) contours: SMA continuum emission at 1.3 mm towards the IRAS 20343+4129 
region, obtained with natural weighting. Contours are -3, 3, 6, 9, 12, 15, and 18 times the rms noise, 2 mJy beam -1 . The synthesized 
beam, shown in the bottom right corner, is 3'.'5 X 2'.'6, at P.A. = -38?0. Thin grey contours: 3.6 cm emission obtained with the VLA 
(Sridharan et al. 2002). Contours are 3, 6, and 9 times the rms noise, 0.2 mJy beam -1 . Grey scale: H2 emission (continuum plus line) 
at 2.12 /urn from Kumar et al. (2002). The crosses correspond to infrared sources from the 2MASS Point Source Catalog (PSC). 



obtain a dust temperature for the cold component of 44 K. 
However, in this last estimate of Tj, the flux densities were mea- 
sured with single-dish telescopes with angular resolutions be- 
tween 10" and 100", including the contribution from IRS 1, 
IRS 3 and other sources in the region. Thus, r<j = 44 K can 
be considered an upper limit. We adopt the intermediate value of 
T d = 30 K. 

The mass of MM7 is -0.7 Mq, and the total mass of MM1 to 
MM4 is ~ 6 Mq. As a reference, the masses derived for the east- 
ern and western single-dish condensations are 44 and 23 Mq, re- 
spectively (Beuther et al. 2002a, 2005; the authors adopt a dust 
temperature and an opacity law that yield masses very similar 
to those obtained using our assumptions). Thus, the millimeter 
compact sources detected with the SMA are embedded in a more 
massive gas halo. For MM6, the dust condensation associated 
with IRS 1, we obtained a mass of ~ 0.8 Mq. There is one con- 
densation, MM5, 5" to the north of IRS 1 that is slightly offset 
(~ 2") to the west of the infrared source IRS IN (Fig.Q]). Since 
the SMA and 2MASS positional uncertainties are 075 and 076 
(Skrutskie et al. 2006), respectively, it is not clear from our data 
whether MM5 is associated with IRS IN or is tracing a different 
source. Finally, we did not detect IRS 3 at 1.3 mm, setting an 
upper limit for its cirmcumstellar mass of ~ 0.2 Mq. It is worth 



noting that IRS 3 is the only source associated with centimeter 
continnum emission in the field (see Fig.Q] and § 14.21 ). 

3.2. CO (2-1) 

Channel maps of the CO (2-1) emission are displayed in Fig. [2] 
CO (2-1) emission extends from -2 up to 33 km s _1 , with the 
systemic velocity being 11.5 km s _1 . The strongest CO (2-1) 
feature is associated with IRS 1, and spans several channels for 
blueshifted and redshifted velocities, as can be seen also in the 
spectrum of the CO (2-1) emission toward IRS 1 (Fig.|2]i. 

The map of the low-velocity emission, integrated from 8.4 
to 14.8 km s _1 , is shown in Fig. [3] From the figure, we can see 
clearly the association of CO (2-1) with IRS 1, as well as the 
presence of low-velocity components southwards of IRS 1 and 
surrounding IRS 3 that are associated with the large-scale CO(2- 
1) blueshifted lobe from Beuther et al. (2002b). It is worth noting 
that the two elongated structures on either side of IRS 3 seem 
to be associated with the fan-shaped structure found in H2 by 
Kumar et al. (2002). 

Regarding the high velocities, these are only present in the 
immediate surroundings of IRS 1. Figure [4^ plots the integrated 
high-velocity emission toward IRS 1 . Blue velocities have been 
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Fig. 2. CO(2-l) channel maps of the IRAS 20343+4129 region, averaged over 2.11 km s _1 wide velocity intervals. The central 
velocity of each channel is indicated in the upper left corner, and the systemic velocity is 11.5 km s _1 . The crosses indicate the 
position of IRS 1 (north) and IRS 3 (south). The synthesized beam, shown in the lower right corner, is 3'.'5 x 2'.'7, at P.A.= -37?8. 
Contours are -6, -3, 3, 6, 9, 15, 30, and 40 times the rms noise, 0.4 Jy beam -1 . The lower right panel is the spectrum of the CO (2-1) 
emission toward the position of IRS 1. The conversion factor is 2.46 K (Jy/bearrfT 1 (in the Rayleigh- Jeans assumption). 



integrated from -6.4 to 8.5 km s , and red velocities from 14.7 
to 32.8 kms -1 . The high-velocity CO (2-1) emission has abipo- 
lar structure, with the center at the position of IRS 1 , and is elon- 
gated in the east-west direction. Note that the red lobe splits up 
into two subcomponents. 



A position-velocity (p-v) plot obtained toward IRS 1 in the 
east- west direction is shown in Fig.[4j5. Up to ±8" from the zero 
offset position, the emission shows a bipolar morphology, reach- 
ing high-velocities that are blueshifted for positive offsets (to 
the east), and redshifted for negative offsets (to the west). The 
distance from IRS 1 where we find high velocities allows us 
to constrain whether such velocities are due to gravitationally 
bound motions. We find weak high-velocity gas at ~ 12 km s~' 
offset from the systemic velocity at position offsets up to 8", 
or 10000 AU. Such velocities at these distances would imply 
an extremely high central mass of ~ 1000 Mq for the motions 



to be gravitationally bound. Hence the bipolar structure seen in 
CO (2-1) toward IRS 1 is most likely tracing an outflow motion. 

Additionally, we computed the p-v plot toward IRS 1 in the 
north-south direction (Fig. |4j;). The only clear CO feature along 
the direction of the cut is the clump at the offset position zero, 
and does not seem related to any other CO feature to the south 
of IRS 1 , although one would expect a north-south bipolar struc- 
ture judging from the single-dish CO map from Beuther et al. 
(2002b). The CO emission at position zero arises from IRS 1, 
with the high-velocity component coming from outflow motions 
(see above). Note that observing outflow emission in both the 
east-west (Fig.UJi) and north-south (Fig.|4j;) directions could be 
indicating that the IRS 1 outflow is a wide-angle outflow (this 
could be also an effect of not resolving the base of the out- 
flow, but the high-velocity emission of Fig. |4j;, specially in the 
blueshifted lobe, is partially resolved). Finally, the decrease in 
CO emission at around the systemic velocity near the zero off- 
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set position, also seen in the CO spectrum of Fig. [2] may be 
produced by a combination of the missing short spacings and 
opacity effects. Given that the brightness temperature at the line 
peak is around 20 K, similar to the kinetic temperature, the line is 
optically thick at systemic velocities. The CO emission at these 
velocities is probably self-absorbed by foreground quiescent ma- 
terial of the cloud in which the YSOs are embedded. 

We calculated the energetics of the outflow associated with 
IRS 1 for each lobe separately, and listed the values in Table [2] 
The expression used for calculating the outflow CO column den- 
sity jV(CO) from the transition J —> J — 1 (derived from Eq. Al 
of Scoville et al. 1986) is: 



+41°40'30 



N(CO) 



4.33 x 10 1 



r- 



exp 



2.777(7+ 1) 



1 - e- T « 



f T B (v)dv 



K km s 



where T ex is the excitation temperature, to is the optical depth, 
and Tb(v) is the brightness temperature profile. 

For the mass derived from CO, we adopted a mean molecular 
weight per H2 molecule of 2.8, and a CO abundance X(CO) = 
1Q- 4 (Scoville et al. 1986): 



M 

Mo 



2.25 x 10 



-16 



pc z 



JV(CO) 



with A being the area of the line emission. 

We assumed optically thin emission in the line wing, and an 
excitation temperature of ~ 25 K, estimated from the spectrum 
in Fig. 12 Due to the lack of observations of other CO transitions, 
we could not make a better estimate of T ex . However, this effect 
is small, as varying T ex between 15 and 30 K yields to a vari- 
ation in the outflow parameters of only ~ 7%. For the red lobe 
we integrated from 15 to 33 km s -1 , and for the blue lobe from 
-6 to 8 km s -1 . The age or dynamical timescale t^ yn was derived 
by dividing the size of each lobe (from the first contour shown 
in Fig. 2k) by the maximum velocity reached in the outflow with 
respect to the systemic velocity (21.5 km s~' for the red lobe, 
and 17.5 km s -1 for the blue lobe). We did not correct for the 
inclination angle, since this parameter is not well known. To ap- 
ply this correction, the velocity must be divided by sin i, and the 
linear size of the lobes must be divided by cos i, with i being the 
inclination angle with respect to the plane of the sky. 

3.3. Infrared emission from 2MASS 

We extracted a sample of infrared sources within the SMA pri- 
mary beam toward IRAS 20343+4129 from the 2MASS Point 
Source Catalog (PSC, Skrutskie et al. 2006), with the aim of 
finding the possible infrared sources associated with the cloud 
of gas and dust studied in this work, and plotted their (J - H), 
(H-K) diagram (Fig.[5k). In the diagram, there are three sources 
with low values of the color indices, including IRS 3. These are 
unreddened (or only slightly reddened) stars. We measured the 
infrared excess as the difference between the (H - K) color and 
the (H-K) color correspoding to a reddened main-sequence star. 
There is a group of five sources, listed in Table [3] for which the 
infrared excess is larger than 1 . We assume that such a large in- 
frared excess (typically, infrared excesses for Class II sources are 
smaller than 0.4; Meyer et al. 1997) is indicative of the YSOs be- 
ing associated with the IRAS 20343+4129 star-forming region. 
Out of these five sources, we detected dust continuum emission 
only toward IRS 1 and possibly IRS IN. 

In order to estimate the spectral type of the infrared sources, 
we plotted them in a J, (J - H) diagram (Fig. ^p). In this di- 
agram, if we deredden IRS 3 along the extinction vector, it 




20 h 36 m 09 s 08 s 07 s 

Right Ascension (J2000) 

Fig. 3. Thick contours: zero-order moment map for the 
low-velocity CO (2-1) emission. Velocities have been inte- 
grated from 8.4 to 14.8 km s _1 . Contours range from 8 to 
56 Jy beam -1 km s" 1 , increasing in steps of 8 Jy beam -1 km s . 
Thin contours: zero-order moment of the blueshifted CO (2- 
1) emission observed in single-dish by Beuther et al. (2002b). 
Contours are 5.2, 6.3, 7.4, 8.5, and 9.6 Jy beam -1 km s _1 . Grey 
scale: H2 emission (continuum plus line) at 2. 12 yum from Kumar 
et al. (2002). Note that the SMA low-velocity structures on ei- 
ther side of IRS 3 are associated with H2 extended emission. The 
crosses correspond to infrared sources from the 2MASS Point 
Source Catalog (PSC). 



falls at the position of B2 stars, consistent with the spectral 
type derived from its centimeter continuum emission (see § 14.21 
and Miralles et al. 1994). Assuming that IRS 3 is a B2 main- 
sequence star, we derived the amount of visual extinction toward 
IRS 3, A v = 13.3 mag. Regarding IRS 4, IRS 5 and IRS 6, they 
all have spectral types around K0 or later. 

In the (J-H), (H-K) diagram of Fig.[5k, when dereddening 
IRS 1 and IRS IN by A v = 13.3 mag (the visual extinction 
toward IRS 3), we found that IRS 1 remains inside the loci of 
luminous Class I and Herbig Ae/Be stars (Lada & Adams 1992; 
Lee et al. 2005), while IRS IN, after dereddening, has colors 
similar to YSOs of low luminosity. 

4. Discussion 

4.1. The young high-velocity bipolar outflow toward IRS 1 

Small scale outflow found toward IRS 1: The parameters of 
the IRS 1 outflow are similar to the mean values of low-mass 
outflows (Wu et al. 2004), and are about 2 orders of magnitude 
smaller than the parameters derived from surveys of high-mass 
outflows (Beuther et al. 2002b; Wu et al. 2005; Zhang et al. 2005; 
Xu et al. 2006). However, given that the surveys of high-mass 
molecular outflows are based on single-dish observations, pick- 
ing up large-scale structure and considering thus larger areas of 
outflow emission, we compared the parameters of the IRS 1 out- 
flow with other high-mass outflows observed with interferome- 
ters, listed in Table [4] From the table we found that the IRS 1 
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Fig. 4. a) CO (2-1) high-velocity emission toward IRAS 20343+4129. Light grey contours, corresponding to the redshifted emission 
in a velocity range from 14.7 to 32.8 km s , range from 2 to 38 Jy beam -1 km s , increasing in steps of 6 Jy beam km s . Black 
contours, corresponding to the blueshifted emission in a velocity range from -6.4 to 8.5 km s _1 , are the same as light grey contours. 
Crosses mark the positions of IRS 1 (north) and IRS 3 (south), b) Position-velocity (p-v) plot along the east-west direction centered 
on IRS 1. Contours are -5, -2, 2, 5, 10, 20, 30, 40, and 45 times 0.4 Jy beam -1 , c) The same as in b) along the north-south direction. 



Table 2. Physical parameters 3 of the outflow driven by IRS 1 





Age 


N u 


Mass 


M 


P 


P 


-^kin 


^mech 


Lobe 


(yr) 


(cm- 2 ) 


(Mq) 


(Mq yr" 1 ) 


(Mq km s" 1 ) 


(Mq km s -1 yr" 1 ) 


(erg) 


(£©) 


Red 


3100 


2.7xl0 16 


0.028 


9.0xl0~ 6 


0.50 


1.6xl0~ 4 


9.0xl0 43 


0.20 


Blue 


3800 


2.6xl0 16 


0.027 


7.2xl0~ 6 


0.38 


l.OxlO- 4 


5.3xl0 43 


0.09 



11 The parameters were obtained as follows. Age: t iyn (see main text); mass-loss rate: M = M/t dyn ; momentum: P = MV mnge (V rangc is the range 
for which we integrated the emission for each lobe, see main text); momentum rate (or mechanical force): P = P/ti yn ; energy of the outflow: 
E = 1/2 M Vl '■; mechanical luminosity: L mcch = E/t iyn . 



outflow is about 2-3 orders of magnitude less energetic than the 
outflows of high-mass protostars observed with high angular res- 
olution. In some cases, the outflow parameters were corrected 
for opacity and inclination, but these effects may contribute only 
about 1 order of magnitude. From this comparison it seems that 
IRS 1 must be a low/intermediate-mass YSO. 

We additionally considered the centimeter continuum lumi- 
nosity produced by shock ionization that should be observed 
for the IRS 1 outflow. From the correlation between the cen- 
timeter luminosity and the outflow momentum rate found by 
Anglada (1995) for a sample of low/intermediate-mass YSOs, 
the IRS 1 outflow can account for a centimeter luminosity of 



~ 0.6 mJy kpc 2 , assuming that all of the stellar wind is shocked. 
This centimeter luminositiy is undetectable with the sensitivity 
of the observations of Sridharan et al. (2002, shown in Fig. [TJ, 
which is on the order of 0.4 mJy kpc 2 . Then, the momentum 
rate derived for the outflow of IRS 1 is not able to produce a 
detectable amount of centimeter continuum emission, which is 
consistent with our observations. 

Finally, the image of the H2 line at 2.12 /im (1-0 S(l)) re- 
veals strong emission very close to IRS 1 (Kumar et al. 2002), 
being elongated in the east-west direction, and thus coincident 
with the direction of the outflow of IRS 1 . This would suggest 
that the H2 emission at 2. 12 jjm close to IRS 1 arises from shocks 




H-K J-H 

Fig. 5. a) (J - H), (H - K) diagram of 2MASS sources within the SMA primary beam toward the IRAS 20343+4129 region. The 
infrared excess is measured with respect to the black thin solid line, which corresponds to the reddening line extending from the 
loci of main-sequence stars, marked with a black thick solid line (Allen 1976, and extinction law of Rieke & Lebofsky 1985). The 
dotted black line indicates the locus of sources with infrared excess equal to one. Note that IRS 3 has only a small infrared excess, 
while IRS 1 and IRS IN are among the sources with highest infrared excess. Thin long-dashed lines mark the region in the diagram 
where most luminous Class I sources and Herbig Ae/Be stars are found (Lada & Adams 1992; Lee et al. 2005). b) J, (J - H) 
diagram of the 2MASS sources in the IRAS 20343+4129 region. The solid black line indicates the locus of the main-sequence stars 
of different spectral types (Allen 1976). Note that when dereddening IRS 3 along the extinction vector (dotted line), it coincides 
with the position of B2 stars, for Ay = 13.3. In both panels, filled circles correspond to the sources with infrared excess larger than 
1. 

Table 3. 2MASS sources within the SMA primary beam with infrared excess > 1 





Identification 


Position 








infrared 


Source 


2MASSJ+ 


a(J2000) 


(5(J2000) 


J 


H 


K d 


excess b 


IRS1 


20360753+4140090 


20:36:07.53 


+41:40:09.1 


15.28 


12.02 


8.90 


1.17 


IRS IN 


20360769+4140121 


20:36:07.69 


+41:40:12.2 


15.18 


13.82 


9.40 


3.61 


IRS3 C 


20360725+4139528 


20:36:07.25 


+41:39:52.8 


12.37 


11.06 


10.13 


0.15 


IRS4 


20360656+4140167 


20:36:06.57 


+41:40:16.7 


16.56 


15.82 


14.17 


1.21 


IRS5 


20360762+4140024 


20:36:07.63 


+41:40:02.5 


15.43 


15.12 


11.79 


3.15 


IRS6 


20360769+4139460 


20:36:07.69 


+41:39:46.1 


15.05 


15.17 


13.94 


1.31 



11 The filter is K s , but we write K for simplicity. 

b The infrared excess is measured as the difference between the measured (H - K) color and the {H - K) color correspoding to a reddened 

main-sequence star (Allen 1976, and extinction law of Rieke & Lebofsky 1985). 
c Although IRS 3 does not have an infrared excess larger than 1, we include this source in the table due to its relevance in the paper. 



in the outflow. However, Comeron et al. (2002) detected line 
emission at 2.225 /mi, which could be due to the 1-0 S(0) line 
of H2. This H2 line at 2.225 fim has been found toward Class 
I and flat-spectrum sources (e. g., Doppmann et al. 2005), asso- 
ciated with outflows (Everett et al. 1995; Davis & Smith 1999; 
Caratti o Garatti et al. 2006), and with photon-dissociated re- 
gions (Ramsay et al. 1993; Luhman et al. 1998), but usually the 
H2 line at 2.225 fim is much fainter than the H2 line at 2. 12 fim, 
while this is clearly not the case for IRS 1 (see Fig. 7 of Comeron 
et al. 2002). Ratios of different intensities of H2 lines are used to 
discriminate between shock-excited H2 emission and excitation 
by fluorescence. A detailed analysis of the ratios of high spectral 
resolution observations of different H2 lines would allow us to 
study the mechanism of the H2 excitation, but this is out of the 
scope of this paper. 



Large-scale CO emission: Single-dish observations of CO (2- 
1) toward IRAS 20343+4129 show a blue lobe centered around 
IRS 3, and a red single-peaked lobe around IRS 1 (Beuther et 
al. 2002b). H2 emission in a fan-shaped structure is found to be 
associated with the blue CO large-scale lobe (see Fig. [3]). This 
seems to suggest the existence of an outflow in the north-south 
direction. However, from the SMA data we found no evidence 
of such a north-south outflow (see Fig.|4^, c). 

We compared the SMA CO (2-1) channel maps with the 
single-dish CO (2-1) channel maps of H. Beuther. The blue lobe 
observed with the single-dish data results from integrating only 
from 8 to 9 km s _1 , an interval very close to the systemic ve- 
locity. In the SMA channel maps from Fig. [2] we only detected 
emission at the position of the blue single-dish lobe for veloc- 
ities between 8.4 and 14.8 km s _1 . Considering that the large- 
scale blue lobe arises from extended emission of ~ 25" in size, 
we found that this emission could not have been detected by the 
SMA. Given the shortest baseline of an interferometer, and the 
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Table 4. Physical parameters of the IRS 1 outflow compared with other low, intermediate, and high-mass outflows observed with 
interferometers 







Age 


Mass 


M 


P 


P 






^mech 




Region 


(Lq) 


(yr) 


(M ) 


(M Q yr" 1 ) 


(Mokms- 1 ) 


(M© km s" 1 


yr 1 ) 


(erg) 


(Le) 


Ref. 


HH211 


3.6 


1400 


0.0024 


1.7xl0- 6 


0.040 


2.8x10- 


-5 


6.9xl0 42 


0.027 


1 


121391 


440 




0.14 




3.6 


1.4x10- 


-3 


1.2xl0 45 




2 


I20343-IRS1 


3200 


3400 


0.055 


1.6x10^ 


0.88 


2.6x10" 


-4 


1.4X10 44 


0.29 


3 


I20293-A 


6300 


4300 


2.0 


4.5x10"* 


90 


2.1x10- 


-2 


4.1xl0 4b 


79 


4 


118182 


20000 




7.3 




71 






8xl0 45 




5 


ON2N 




37000 


58 




1060 


2.8x10- 


-3 


2.0xl0 47 


45 


6 



REFERENCES: 1: Palau et al. (2006); 2: Beltran et al. (2002); 3: this work; 4: Beuther et al. (2004a); 5: Beuther et al. (2006); 6: Shepherd et 
al. (1997). 



size of the observed emission, one can estimate the fraction of 
correlated flux detected by the interferometer. For a source of 
25", this corresponds to a 0.2% for our SMA configuration (with 
a shortest baseline of 1 1 kA). Thus, out of the total CO (2-1) flux 
observed in single-dish for the blue lobe, 1310 Jy km s _1 , the 
SMA could pick up only 2.6 Jy km s _1 , which should be de- 
tected at a few times the rms of the SMA channel maps, as is 
the case (see channel at 9.5 km s _1 of Fig. |2j. We concluded that 
almost all of the large-scale blueshifted lobe seen in single-dish 
has been resolved out by the SMA. Note that in the single-dish 
images of Beuther et al. (2002b) there is a blue contour at the po- 
sition of IRS 1, suggesting that if blueshifted emission had been 
integrated for velocities < 8 km s , the blue lobe of the IRS 1 
outflow detected with the SMA would appear in the single-dish 
observations as well. 

Regarding the single-dish red lobe, this was obtained inte- 
grating from 13 to 15 km s _1 (Beuther et al. 2002b). We esti- 
mated the contribution of the SMA redshifted emission to the 
single-dish redshifted emission. The SMA integrated intensity 
of the red lobe for the same range of velocities as Beuther et al. 
(2002b) is 160 Jy km s , obtained from the spectrum toward 
IRS 1, and assuming a size of the lobe of ~ 10". Beuther et 
al. (2002b) obtain a flux density of 1300 Jy km s , implying 
that the SMA redshifted emission can account for about 12% 
of the redshifted emission detected with the single-dish data, a 
contribution about 2 orders of magnitude higher than that for the 
blueshifted large-scale lobe. This suggests that the blue and red 
lobes seen in the single-dish data, come from material at differ- 
ent spatial scales. 

4.2. Is IRS 3 driving a cavity around it? 

As seen above, the CO (2-1) emission from the single-dish data 
(Beuther et al. 2002b) shows a slightly blueshifted large-scale 
lobe that corresponds well with the fan-shaped structure seen in 
H2 emission around IRS 3 (see Fig. [3). In addition, from the 
SMA data of this work, we found low-velocity CO elongated 
structures on either side of IRS 3, as well as dust condensations 
also on either side of IRS 3. All these observational features 
seem to suggest that IRS 3 is interacting with the surrounding 
medium and producing a shell of circumstellar gas expanding 
away from IRS 3. We search for any kinematical evidence of 
such an expanding cavity in the low-velocity CO emission, but 
the strong missing flux, and the presence of at least one bipo- 
lar outflow and other YSOs surrounding IRS 3, make this search 
difficult. However, a preliminary study of the NH3 emission in 
the region revealed dense gas surrounding IRS 3 whose kine- 



Table 5. Summary of the centimeter observations carried out to- 
ward IRAS 20343+4129 



A 




Observing 


#FWHM 




(cm) 


(mJy) 


year 


(") 


Refs. 


6 


1.1+0.2 


1989 


5" 


Miralles et al. (1994) 


3.6 


1.3 + 0.3 


1994 


8" 


Carraletal. (1999) 




1.8 + 0.1 


1998 


1" 


Sridharan et al. (2002) 


2 


< 1.0" 


1989 


5" 


Miralles et al. (1994) 


0.7 


< 9 a 


2003 


1'.'5 


Menten et al., in prep. 



a Upper limits at the 3<x level. 



matics were consistent with an expanding shell with IRS 3 in 
its center (Palau et al., in prep.). In the following, we consider 
whether IRS 3 is able to drive such a cavity. 

Given that the centimeter continuum emission associated 
with IRS 3 (of 1.8 mJy, Sridharan et al. 2002) is consistent with 
an ionizing B2 star (Panagia 1973), and that the luminosity-to- 
mass ratio of typical B2 stars is high enough to allow the star to 
push away the surrounding material by radiation pressure (given 
typical dust opacities of molecular clouds, Calvet et al. 1991; 
Anglada et al. 1995), the first interpretation to explore is that the 
cavity is driven by the radiation pressure from IRS 3. This leads 
to a scenario in which IRS 1 is a low/intermediate-mass YSO 
accounting for a small fraction of the total bolometric luminos- 
ity of IRAS, while IRS 3 is a high-mass star accounting for most 
of the bolometric luminosity. However, the cavity could also be 
driven by a stellar wind from IRS 3. We studied the scenario of 
a wind-driven cavity by following the model of Anglada et al. 
(1995). In this model, the centimeter continuum emission is as- 
sumed to trace an ionized stellar wind. We adopted a spectral 
index for the wind of 0.6, and estimated a mass loss rate of the 
ionized material following Beltran et al. (2001). By noting from 
the observations the radius of the cavity (~ 10"), the velocity of 
its walls (~ 2 km s" 1 ), and the external pressure of the ambient 
cloud (~ 3.3 km s , from the line width of the CO (2-1) line 
in the walls of the cavity), we obtained a maximum radius, den- 
sity, and dynamical timescale of the cavity of 12", 2300 cm -3 , 
and 11000 yr, respectively, which are reasonable values for the 
region. This leads to a scenario in which IRS 3 is not necessar- 
ily the most massive source of the region, and then IRS 1 could 
account for most of the IRAS luminosity. 

Thus, in the interpretation of the centimeter continuum emis- 
sion as either an UCH 11 region, or an ionized stellar wind, IRS 3 
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can create a cavity of swept up material around it. In order to 
distinguish between these two scenarios, it would be very use- 
ful to consider the spectral energy distribution in the centimeter 
range for IRS 3. In addition to the 3.6 cm measurement from 
Sridharan et al. (2002), shown in this work, there are other ob- 
servations at 6, 3.6, and 2 cm, listed in Table [5] The only simul- 
taneous measurements are those at 6 and 2 cm, which result in 
a spectral index < -0.1, consistent with an optically thin UCHn 
region. However, the other observations at 3.6 cm are not consis- 
tent with this flat spectral index. With the measurements at 3.6 
and 6 cm, the spectral index ranges from 0.3 to 0.9 (depending 
on the epoch of the observations at 3.6 cm). This spectral index 
is consistent with emission from an ionized wind (e. g., Panagia 
& Felli 1975). Thus, from the current data it is not clear whether 
the source is variable with time, and what is the value of the 
spectral index of the centimeter source associated with IRS 3. 
New observations at 6, 3.6, 2, and 1.3 cm toward IRS 3 would 
give insight into which mechanism, a stellar wind or radiation 
pressure, is driving the cavity. 

4.3. On the nature of IRS 1 

From an analysis of JHK color-magnitude diagrams including 
the high-mass stars of the Cygnus OB2 association, Comeron 
et al. (2002) propose that IRS 1 may be one of the most lu- 
minous and embedded objects in the entire association. In the 
K, (H - K) diagram, IRS 1 is the reddest object of the asso- 
ciation, and dereddening along the extinction vector yields a 
very bright K magnitude compared with the other massive stars. 
However, this plot sets an upper limit to the intrinsic brightness 
because in the K band there may be some contribution from 
circumstellar material. For this reason, Comeron et al. (2002) 
use the H, (J - H) diagram, which is not seriously affected 
by circumstellar emission, and find again that IRS 1 is among 
the brightest. Since the color-magnitude diagrams were made 
by using the Second Incremental Release of the 2MASS Point 
Source Catalog (PSC), and this release is by now obsolete, we 
redid the diagrams with magnitudes from the current release of 
2MASS PSC, and found the same values except for the J magni- 
tude, yielding (J-H) = 3.26 instead of 4.23 used by Comeron et 
al. (2002). Thus, dereddening along the extinction vector in the 
H, (J-H) diagram does not set IRS 1 among the brightest mem- 
bers of the association, but yields that IRS 1 must be still intrin- 
sically brighter than stars with spectral type B0, for which typ- 
ical luminosities are around 25000 Lq (Panagia 1973), at least 
one order of magnitude higher than the bolometric luminosity of 
IRAS 20343+4129. 

If IRS 1 were a high-mass ZAMS star, one would expect 
to detect centimeter continuum emission from ionized gas, and 
to detect the star at optical wavelengths. However, no centime- 
ter continuum emission was detected above the 3<x level of 
0.6 mJy beam 1 , and IRS 1 does not appear in the POSSII plates. 
If IRS 1 were a high-mass protostar, deeply embedded and still 
accreting most of its mass, one would expect to observe a mas- 
sive envelope, with a mass of the order of the accreted mass, 
surrounding the protostar. But we did not detect a massive en- 
velope toward IRS 1 from the millimeter continuum emission 
(the circumstellar mass was ~ 0.8 Mq, § 13.11 ), and the outflow 
parameters were comparable to those of low/intermediate-mass 
protostars (§ 14. 11 1. Thus, a high-mass nature for IRS 1 is not con- 
sistent with our observations. All this suggests that IRS 1 cannot 
be considered a reddened stellar photosphere but a YSO with a 
cold envelope, and that the estimation of its brightness from the 
magnitude-color diagrams gives only an upper limit. Note how- 



ever that the JHK magnitudes of IRS 1 cannot be accounted for 
by a low-mass YSO at the distance of the region. 

At this point we consider whether the continuum emission 
is really tracing all the dust surrounding IRS 1 . We considered 
four different possibilities. First, the SMA could be filtering out 
large-scale emission; we considered this option and we ruled 
it out because the largest angular scale (FWHM) to which the 
SMA is sensitive is ~ 8", larger than ~ 4", the observed size 
of the millimeter source associated with IRS 1 (MM6). In ad- 
dition, the CO emission observed with the SMA toward IRS 1 
shows that the SMA is really sensitive to scales larger than the 
size of the dust condensation associated with IRS 1, and that the 
dust emission is significantly more compact than the CO emis- 
sion. Second, the dust could be optically thick; however, this 
possibility yields an unrealistically low value for the dust tem- 
perature, given the flux density at 1.3 mm and the deconvolved 
size for MM6 of 3". Third, the dust could be sublimated; since 
dust sublimates at ~ 1500 K (Whitney et al. 2004), a consider- 
able amount of gas at such a high temperature should be detected 
at optical wavelengths. Fourth, CO emission has been detected 
with no continuum emission for a few YSOs; however, these ob- 
jects are evolved (Class II/III) low-mass systems (e. g., Andrews 
& Williams 2005; Takeuchi & Lin 2005). Thus, none of these 
possibilities is convincing for the case of IRS 1, hinting that the 
mass traced by the millimeter continuum emission is most likely 
all the circumstellar mass associated with IRS 1 . 

Therefore, the current data suggest that IRS 1 is a 
low/intermediate-mass YSO. In order to further constrain the 
mass and the evolutionary stage of IRS 1, we plotted the spectral 
energy distribution (SED) compiled from 2MASS (corrected for 
interstellar extinction, see § 13.31 ). MSX, IRAS, observations in 
the submillimeter range carried out with SCUBA on the JCMT 
(Williams et al. 2004), and the SMA (this work). As the IRAS 
flux densities may have contribution from both IRS 1 and IRS 3, 
they are only upper limits. The fluxes measured by SCUBA are 
also upper limits because the single dish is picking up large-scale 
emission, partially arising from IRS 3 and from the dust conden- 
sations on either side of IRS 3. The angular resolution of the 
MSX images allowed us to estimate the flux density from IRS 1 
by integrating the mid-infrared emission in an aperture of ~ 15" 
of diameter around IRS 1 . 

The resulting SED (Fig. |6]) shows a steep profile for the 
2MASS wavelengths, suggesting that any contribution from a 
hot photosphere is negligible. Rather, the peak of the SED lying 
between 10 and 100 fim indicates that most likely a cold enve- 
lope dominates the SED, with d(log/lf / i)/d(log/l)> between 
2 and 100 /vm, consistent with the classification of IRS 1 as a 
Class I source (e. g., Hartmann 1998). Note that if IRS 1 were a 
Class source, the SED should not show significant emission in 
the near infrared (Andre et al. 1993; Lada 1999). We compared 
the IRS 1 SED with SEDs of Class I sources of bolometric lumi- 
nosities between 4 and 960 Lq from the literature, scaled to the 
distance of IRAS 20343+4129. From this comparison, we found 
that IRS 1 is likely not a low-mass source, but rather its SED 
resembles that of intermediate-mass sources with bolometric lu- 
minosities around 1000 Lq (the exact value of the luminosity 
of IRS 1 cannot be determined because the IRAS flux densities 
are only upper limits, and could be the main contribution to the 
bolometric luminosity). From the evolutionary tracks of Palla & 
Stahler (1993) for intermediate-mass stars, and assuming that 
IRS 1 is already in the birthline, one would expect a stellar mass 
for IRS 1 of around 5 Mq. This is a factor of 6 higher than the 
circumstellar mass derived from the SMA data, ~0.8 Mq (§ 13. U . 
In the low-mass regime, Class sources are expected to have 
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Fig. 6. Spectral Energy Distribution (SED) for IRS 1. Black dots cor- 
respond to IRS 1, and thin curves are SEDs from the literature for 
Class I sources of 960 Lq (dotted, VMR-D-IRS 13: Massi et al. 1999), 
250 L Q (solid, VMR-D-IRS 14: Massi et al. 1999), and 4 Lq (dashed, 
IRAS 04016+2610: Eisner et al. 2005), which have been scaled to the 
distance of IRAS 20343+4129. IRAS and SCUBA fluxes are upper lim- 
its because it is not possible to disentangle the contribution of IRS 1 
from the surrounding sources. For comparison, open triangles corre- 
spond to IRS 3 2MASS fluxes. 2MASS fluxes have been corrected for 
interstellar extinction (see § 13. 3b . 



circumstellar/envelope masses similar to the stellar masses, and 
thus one could argue that 0.8 Mq is too low for an intermediate- 
mass YSO of 1000 Lq. However, one expects that in the Class I 
stage, the ratio of the circumstellar-to-stellar mass progressively 
decreases. Taking into account that intermediate-mass sources 
evolve faster to the main-sequence than the low-mass sources, it 
may be reasonable to have a circumstellar mass lower than its 
stellar mass. In addition, there are some cases in which the mass 
of the dust/gas condensation associated with intermediate-mass 
YSOs is significantly lower (at least by 1 order of magnitude) 
than the estimated stellar mass (Beuther et al. 2004b; Martin- 
Pintado et al. 2005; Zapata et al. 2006). Thus, a circumstellar 
mass of ~ 0.8 Mq for IRS 1 seems to be compatible with its 
classification as a Class I source of around 1000 Lq. 

4.4. Sources in different evolutionary stages: millimeter vs 
near-infrared emission 

In order to gain insight into the star formation process in the 
region, we tentatively classified the different sources identified 
in this work in three different evolutionary stages, depending on 
their millimeter and infrared emission. 



Sources at the end of the accretion phase: We classified in 
this evolutionary stage IRS 3 to IRS 6 (Table [3]), the sources 
detected only in the infrared. From the centimeter and infrared 
emission, IRS 3 seems to be a B2 star. As for IRS 4 to IRS 6, they 
all have spectral types around K0 or later (as shown in Fig.[5t>), 
and thus are low-mass YSOs. Since we did not detect emission 
toward these sources at 1.3 mm, we estimated an upper limit for 
their associated mass of ~ 0.2 Mq (with the same assumptions 
as § 13. 11 1. These YSOs are possibly Class II/III sources. 



Sources in the main accretion phase: We included in this 
group the sources showing both infrared and millimeter emis- 
sion, as is the case of IRS 1. In § !3.3l and l4.3l we concluded that 
IRS 1 seems to be an intermediate-mass Class I source. Another 
source that could be in the main accretion phase is IRS IN, if we 
assume that it is associated with MM5 (see § 13. 11 1. For this source 
we estimated a spectral type AO or later, from the magnitude- 
color diagram of Fig. \5&, and an associated mass of ~ 0.9 Mq, 
from the SMA data. Thus, if IRS IN is an embedded YSO, it 
is a low-mass object. However, IRS IN is the source with the 
highest H - K color in the region (see Fig.|5^), and a possibility 
for such a high H - K color is that the /f-band is contaminated 
by the H2 line at 2.12 /mi, as the continuum-subtracted H2 im- 
ages from Kumar et al. (2002) suggest. Thus, it remains unclear 
from the available data whether IRS IN is a low-mass embedded 
YSO, or an infrared source whose emission mainly arises in the 
interaction of an outflow (either the IRS 1 outflow, or an outflow 
from IRS 1N/MM5 itself) with the surrounding medium. 

Starless core candidates: The sources MM1 to MM4 and 
MM7 (Table [TJ have been detected only in the millimeter, and 
lie in a region completely dark in the near infrared (Fig. [TJ. The 
masses of these dust condensations are 0.7-1.2 Mq, and thus 
they are low-mass condensations. Sources bright only in the 
millimeter could also be tracing Class protostars being at the 
beginning/main accretion phase. However, given that we did 
not find any sign of star formation such as outflow emission, we 
suggest that some, if not all, of these sources could be starless 
core candidates. 

Therefore, the IRAS 20343+4129 region harbors sources 
that seem to be in different evolutionary stages, and with stellar 
masses ranging from < 0.2 to 8-10 Mq. The intermediate/high- 
mass sources of the region, IRS 1 and IRS 3, are in different 
evolutionary stages. IRS 3 is visible at optical wavelengths, does 
not show infrared excess, and has no CO (2-1) nor dust emis- 
sion associated, suggesting that it has already finished the main 
accretion phase. On the contrary, IRS 1 is not detected in the 
optical, has strong infrared excess and dust emission associated, 
and is the driving source of a bipolar outflow, which indicates 
that IRS 1 is still accreting a significant amount of mass, with 
IRS 1 therefore younger than IRS 3. However, even with dif- 
ferent evolutionary stages, IRS 1 and IRS 3 could have formed 
simultaneously, since the uncertainty in their masses is impor- 
tant, and the contraction time to the main-sequence (since the 
formation of a first hydrostatic core until the hydrogen burn- 
ing stage), which deacreases with the mass of the YSO, in the 
intermediate/high-mass regime gets comparable or shorter than 
the free-fall timescale, and thus an intermediate-mass Class I 
YSO may have formed simultaneously with a high-mass ob- 
ject already in the main-sequence phase. For example, using 
the computations of Bernasconi & Maeder (1996), and assum- 
ing that IRS 3 is around 9 Mq and IRS 1 is around 5 Mq, 
IRS 3 has a contraction time of -0.3 Myr, and IRS 1 of ~ 1 Myr. 
This is different from the low-mass regime, where YSOs evolve 
to the main-sequence in a timescale > 10 Myr (e.g., Hayashi 
1961; Iben 1965), which is always much larger than the free- 
fall timescale (~ 0.5 Myr). Thus, a low-mass YSO in the pre- 
main-sequence phase (Class II/III) and a low-mass YSO in the 
main accretion phase (Class 0/1) cannot have formed simultane- 
ously. As seen in this section, we find low-mass YSOs in the 
very first phases of formation, and others that seem to be Class 
II/III candidates, and thus they cannot be coeval. Therefore, the 
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star formation in IRAS 20343+4129 seems to be a continuous 
process. It is worth noting that the intermediate/high-mass star 
IRS 3 is in an evolutionary stage equal or more evolved than 
the low-mass sources. This is in contrast to previous claims that 
intermediate/high-mass stars in clusters appear less evolved than 
the low-mass stars in the cluster (Massi et al. 2000; Kumar et al. 
2006). 

Finally, the spatial distribution of the sources found in 
IRAS 20343+4129 shows that star formation in this region is 
localized, that is, to the north and south of IRS 3 there are YSOs 
that are already bright in the infrared, while to the east and to the 
west there are millimeter sources that seem to be starless cores. 
Since IRS 3 and IRS 1 could be coeval, we suggest that these 
two sources are reflecting the initial conditions of high density 
in the parental cloud. As for the low-mass dust condensations on 
either side of IRS 3, we propose that they may be compressed by 
the expanding cavity driven by IRS 3, and thus could eventually 
form a new generation of stars. 
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Fig. 7. Grey scale: H2 emission (continuum plus line) at 2.12 
fan in the IRAS 20343+4129 region (Kumar et al. 2002). 
White contours correspond to the 3.6 cm emission from 
Carral et al. (1999), and are 0.66, 0.99, 1.14, 1.33, 1.52, and 
1.71 mJy beam -1 . Black contours are the 1.2 mm continuum 
emission observed in single dish by Beuther et al. (2002a), 
and trace the dust cloud where IRS 1 and IRS 3 are forming. 
Contours are: 25, 50, 100, 200, and 300 mJy beam -1 . Note that 
the centimeter emission traces the ionized southwestern border 
of the cloud, as well as IRS 3 in the center of the cloud. Crosses 
indicate the spatial distribution of the 2MASS sources in the re- 
gion with infrared excess larger than zero: (0-0.4) for the yellow 
color (Class II sources); (0.4-1) for the orange color (Class I), 
and > 1 for the red color, following Matsuyanagi et al. (2006). 
Note that most of the 2MASS sources with infrared excess > 1 
lie in the center of the cloud, near IRS 3 and IRS 1 . 



4.5. IRAS 20343+4129 within Cygnus OB2 

The H2 emission at 2.12 /jm shows a cometary arch about 1' 
to the southwest of IRS 3, which is also detected in centime- 
ter emission (see Fig.|7j Kumar et al. 2002; Carral et al. 1999). 
This cometary arch follows the border of a dust cloud traced 
by 1.2 mm emission observed in single-dish by Beuther et al. 
(2002a), and could be produced by the ionization front from 
a nearby OB star. In fact, the arch is facing the center of the 
Cygnus OB2 association. This would be similar to the bright- 
rimmed clouds facing HII regions, such as IC 1396N (Sugitani 
et al. 1991; Beltran et al. 2002), which is ionized by an 06.5 
star at 13 pc of distance (Schwartz et al. 1991). We plotted the 
OB stars of the association (Reed 20020), and found that there 
are five 06-09 stars about 25' (9 pc at the adopted distance 
for IRAS 20343+4129) to the south-west of the arch. Given 
the flux density of the centimeter continuum emission tracing 
the cometary arch, one can estimate the required flux of ion- 
izing photons per unit area that must reach the cloud to pro- 
duce the observed centimeter emission (Lefloch et al. 1997). 
In our case, the flux density of the arch at 3.6 cm is around 
5 mJy, and this requires a flux of ionizing photons per area unit 
of 4 x 10 9 s -1 cirT 2 . The flux of ionizing photons for an 06 
star is tipically 1.20 x 10 49 s -1 (Panagia 1973), and assuming 
a distance to the arch of 9 pc, the flux of ionizing photons per 
area unit reaching the cloud is 1.2 x 10 9 s -1 cm -2 , close to the 
value required to account for the flux of the centimeter emis- 
sion in the arch. Thus, the 06 star 25' to the southwest of the 
arch, BD+41 3807, is probably the star ionizing the border of 
the IRAS 20343+4129 cloud. Finally, we calculated the infrared 
excess (as described in § |3.31 > of the 2MASS sources within a 
diameter of 4' centered on IRS 3 (thus including a region out- 
side the rim), and plotted the sources with positive infrared ex- 
cess in Fig. [7] In the figure, crosses with redder colors corre- 
spond to sources with larger infrared excess (and thus presum- 
ably younger). As seen from the figure, most of the sources with 
strongest infrared excess within the field are found inside the 
cloud, thus constituting a region of recent star formation as com- 
pared to its surroundings. In addition, the millimeter sources de- 
tected in this work could be sites of future star formation, and 
thus IRAS 20343+4129 is a region actively forming stars within 
the Cygnus OB2 association. 



5. Conclusions 

We observed the dust continuum emission at 1.3 mm with the 
SMA as well as the CO (2-1) emission toward the massive star- 
forming region IRAS 20343+4129, in order to study the prop- 
erties of the different protostars in the region and their interac- 
tions with the surrounding medium. Two bright infrared sources, 
IRS 1 in the north and IRS 3 in the south, lie inside the SMA pri- 
mary beam, and IRS 3 is associated with centimeter continuum 
emission. Our main conclusions can be summarized as follows: 

1 . The dust continuum emission reveals three main condensa- 
tions, to the north (associated with IRS 1), to the east and 
to the west of IRS 3, with the western condensation being 
the brightest one and consisting of different subcondensa- 
tions, of ~ 1 Mq each. Toward the eastern condensation, 
of ~ 0.7 Mq, the SMA has filtered out most of the emis- 
sion. The estimated mass of the condensation associated with 
IRS lis -0.8 M Q . 



3 Catalog available at http://othello.alma.edu/~reed/OBfiles.doc 
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2. We discovered a bipolar high-velocity CO outflow, elongated 
in the east-west direction, and identified IRS 1 as its driving 
source. The millimeter continuum and CO emissions indi- 
cate that IRS 1 is not a high-mass YSO, and the SED agrees 
with IRS 1 being an intermediate-mass Class I source. 

3. The low-velocity CO emission shows, in addition to emis- 
sion toward IRS 1, two elongated structures on either side of 
IRS 3, coincident with extended H 2 emission at 2. 12yum. The 
emission from the blueshifted lobe of the large-scale CO out- 
flow seen in single-dish by Beuther et al. (2002b) has been 
filtered out by the SMA, and the lobe does not have com- 
pact emission. A scenario in which the blue large-scale CO 
lobe traces a cavity blown up by IRS 3, and where the dust 
condensations on either side of IRS 3 are the result of the 
accumulation of mass in the walls of the expanding cavity, is 
consistent with the observations. In this scenario, the elon- 
gated low-velocity CO emission and the H2 extended emis- 
sion trace the walls of the cavity. The expanding cavity could 
be either driven by a stellar wind from IRS 3, or driven by 
radiation if we assume that IRS 3 is a B2 star. 

4. We found objects with very different properties and evolu- 
tionary stages that have been born in the same parental cloud. 
In addition, these objects are not randomly distributed in the 
cloud, but their distribution seems to be determined, at least 
partially, by the accumulation of mass at the walls of an ex- 
panding cavity driven by IRS 3. 
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